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In living cells, membrane morphology is regulated by various proteins. Many membrane reshaping 
proteins contain a Bin/Amphiphysin/Rvs (BAR) domain, which consists of a banana-shaped rod. 

The BAR domain bends the biomembrane along the rod axis and the features of this anisotropic 
bending have recently been studied. Here, we report on the role of the BAR protein rods in inducing 
membrane tubulation, using large-scale coarse-grained simulations. We reveal that a small sponta¬ 
neous side curvature perpendicular to the rod can drastically alter the tubulation dynamics at high 
protein density, whereas no significant difference is obtained at low density. A percolated network is 
intermediately formed depending on the side curvature. This network suppresses tubule protrusion, 
leading to the slow formation of fewer tubules. Thus, the side curvature, which is generated by 
protein-protein and membrane-protein interactions, plays a significant role in tubulation dynamics. 

We also find that positive surface tensions and the vesicle membrane curvature can stabilize this 
network structure by suppressing the tubulation. 


The Bin/Amphiphysin/Rvs (BAR) superfamily pro¬ 
teins regulate the membrane shape of cell organella as 
well as membrane fusion and fission; therefore, BAR pro¬ 
tein dysfunction is implicated in neurodegenerative, car¬ 
diovascular, and neoplastic diseases 09- However, the 
manner in which these proteins assemble on the biomem- 
brane and cooperate to reshape the membranes is not 
well understood. The extension of membrane tubes from 
liposomes and specific adsorption of the BAR superfam¬ 
ily proteins onto tube regions have been observed in in 
vitro experiments 0S6- 15^ . Frost et al. have experi¬ 
mentally determined that F-BAR proteins are adsorbed 
on flat regions of lipid membranes using electron mi¬ 
croscopy [9|. Although the assembly seems to constitute 
the nucleus of the tubule formation, the tubule protrusion 
process has not been experimentally observed. Recently, 
Tanaka-Takiguchi et al. reported that the formation dy¬ 
namics of tubules from a liposome can differ significantly 
for different F-BAR proteins 131. That is, FBP17 and 
CIP4 simultaneously generate many tubule protrusions 
over the entire liposome surface, while PSTPIP1 and Pac- 
sin2 generate only a few protrusions from a narrow region 
of the surface. In particular, the tubules induced by CIP4 
and PSTPIP1 have the same radius. Thus, the tubule nu- 
cleation process depends on the protein type. However, it 
is not known what causes this difference in tubule nucle- 
ation behaviour. Tanaka-Takiguchi et al. also reported 
that the full length of Pacsin2 induces tubulation, but 
its F-BAR domain region alone does not [l3j]. In con¬ 
trast, Wang et al. reported the tubulation is induced by 
F-BAR domain of Pacsinl more than by the full-length 
protein 0 - 

In the last decade, interactions between laterally 
isotropic objects on biomembrane, such as transmem¬ 
brane proteins and adsorbed spherical colloids, have been 


sives have not yet been explored so far. The BAR do¬ 
mains are banana shaped and generate an anisotropic 
curvature different from the isotropic spontaneous curva¬ 
ture Co 0. This anisotropic nature has recently been 
receiving increasing theoretical interest. The classical 
Canham-Helfrich curvature free energ y Eill . [25l has been 


extended to anisotropic curvatures [26 b 28| . Dommersnes 
and Fournier have derived a many-body potential of long- 
range interactions between point-like anisotropic inclu¬ 
sions and found linear assemblies and egg-carton mem¬ 
brane structures using Monte Carlo simulations 0,0. 
In addition, the adsorption and assembly of BAR do¬ 
mains have been investigated u sing atomic and coarse¬ 
grained molecular simulations |3ll434l |. For example, 
Simunovic et al. have simulated a linear aggregation 
of N-BAR domains parallel to the domain axis 3j, 0. 


However, the relationship between this aggregation and 
tubulation remains unclear. Further, tubular forma¬ 
tion has been simulated using a dynamically triangulated 


membrane model 


and, also, meshless membrane 


intensively investigated }16h22|] . In contrast to such stud¬ 


ies, however, the interactions between anisotropic adhe¬ 


models SHI. Despite these numerous advancements, 
the present understanding of the physics of membrane 
shape deformation due to anisotropic curvature is still 
far from complete. 

In this paper, we focus on the effects of the sponta¬ 
neous (side) curvature C s id e of a protein rod perpendic¬ 
ular to its longest axis on the assembly behaviour. The 
side curvature has not been focused upon in previous 
studies, but here we reveal that it strikingly changes the 
assembly dynamics. The excluded volume or van der 
Waals attraction between proteins and the membrane 
can effectively generate positive or negative C s id e - We 
simulate flat membranes and vesicles using an implicit- 
solvent meshless membrane model 38)-42|, which allows 
a large-scale simulation. A BAR domain is modelled as a 
banana-shaped rod, which is assumed to be strongly ad¬ 
sorbed onto the membrane. The rod length corresponds 
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to r ro d — 20 nm (the BAR domain lengths range from 13 
to 27 nm i)- To investigate the membrane-curvature- 
mediated interactions, no direct attractive interaction 
is considered between the rods. Our previous studies 
showed that parallel and perpendicular assemblies occur 
separately through membrane-mediated attractive inter¬ 
actions at low protein density [38], and that polyhedral 


shapes are formed at high protein density [42J for vesicles 
and membrane tubes. 


Results 

Tubulation from Flat Membrane 

First, we investigate the tubulation from a tensionless 
flat membrane (surface tension 7 = 0 ) at a high rod 
density, (j > IO d = 0.4 (see Fig. [TJ. The protein rods are ini¬ 
tially equilibrated with the rod curvature C ro d = 0 and 
O s ide = 0. Once the spontaneous curvatures are altered 
at t = 0, the rods begin to assemble perpendicularly to 
the rod axis. For a positive spontaneous curvature of 
Osidefrod = 1 (C s id e C ro d > 0 ), many tubules simultane¬ 
ously protrude via the bending of straight rod assemblies 
(see Fig. [fjc) and Supplemental Movie 1). Branches of 
the rod network are formed on the membrane for a short 
time only. When the tubulation is initiated, neighbour¬ 
ing branches are broken through lateral shrinkage of the 
rod assembly. 

For a negative curvature C s idefrod = — 1, the rods form 
a percolated network covering the entire membrane area 
[see the top snapshot in Fig. Ed)] and a tubule pro¬ 
trudes under membrane undulation (see the second snap¬ 
shot in Fig. [ljd) and Supplemental Movie 2). Subse¬ 
quently, tubule growth occurs along the network. Thus, 
the tubulation dynamics is altered remarkably by a rel¬ 
atively small C s ide- Negative and positive C s id e values 
stabilize and destabilize the network branches, respec¬ 
tively. The tubulation at CddeH-od = — 1 is significantly 
slower than that at Cgide^rod = 1 and much fewer tubules 
protrude: the average protrusion time of the first tubule 
are {ttb)/r = 69 ± 4 and 10 ± 1 for CgideA-od = —1 and 1 
at CVod^rod = 4, respectively. 

Such characteristic dynamics is distinguishable from 
the time evolution of the mean cluster size (N c \) and 
the root mean square cluster height ((z as shown 

in Figs. QJe)-(g). For C s id e r r od = -1, the majority of 
the rods belong to one large percolated cluster during 
the tubulation. In contrast, for CsideCrod = 1, {N c \) de¬ 
creases as the tubules are formed and the rod assemblies 
are divided; subsequently, {N c \) slowly increases owing 
to tubule fusion. Branched tubules are formed by this 
fusion [see the bottom snapshot in Fig. QJc) and the late 
stage of Supplemental Movie 1]. Based on the evolution 
of (iV c i), the tubulation pathways are categorized into 
three groups [see Fig. Oik)]: tubulation via percolated- 


network formation (net), via partial-percolated-network 
formation (part), and without percolation (iso). When 
a percolated network does not cover the entire mem¬ 
brane surface or a large cluster of N c \ ~ N ro d is main¬ 
tained for a period shorter than 20 r, we categorize the 
tubulation pathway as part. (A typical dynamics is 
shown in Supplemental Movie 3). For the entire param¬ 
eter range explored in Fig. [ljb), the final structures are 
tubules. As C ro d decreases, the tubulation decelerates 
and a smaller number of large tubules are formed. The 
tubules are nucleated and grow from the network vertices 
at CrodJ'rod = 2.5 or 0.3 and C s ideCod = (see Fig. [2J. 
The tubule radius I?tb is roughly determined by C ro d as 
Rtb ~ 1/Crod- At CVodfrod = 4, the tubule with circum¬ 
ference 27rl?tb — 2 r ro d consists of two hemicylinders of 
the rod assembly. 

Our simulation results show that the network forma¬ 
tion suppresses the tubulation. To confirm this more 
clearly, the effects of C s jde on the rod-membrane interac¬ 
tion are investigated. A percolated network is not formed 
during the tubulation at a low <j > TO d of 0.1. Rather, the 
rods assemble into linear clusters and, subsequently, the 
large clusters (N c \ > 40) transform into tubules (see 
Fig. [3] and Supplemental Movie 4). Although the initial 
cluster formation is slightly slower for negative curvature, 
i.e., Cside^rod = — 1 , no qualitative difference is detected 
in the tubulation dynamics [see Fig.[3][b)]. Thus, we con¬ 
clude that the suppression of the tubulation at high rod 
density is caused by the mesoscale network formation. 

Recent experiments have demonstrated that positive 
surface tension can suppress tubulation by BAR proteins 
f~Hil | and budding by clathrin coats (43|. In our simula¬ 
tion, the positive tension and network formation cooper¬ 
ate to suppress tubulation [see Figs. [U[e)-(g) and[3](b)]. 
At (f> T od = 0.4, the critical tension decreases with increas¬ 
ing C s ide: 7 r rod/^ B F ~ 10 and 70 for C si d e rrod = -1 
and 1 , respectively, where k^T denotes the thermal en¬ 
ergy. These are experimentally measurable magnitudes 
(7 ~ 0.1 mJ/m 2 and 0.7 mJ/rn 2 , respectively). The as¬ 
sembly of rods into a clustered network is not suppressed 
by the applied tension. In contrast, network structure 
breaking does not occur at higher tensions [see Fig.QJf)]. 
Thus, the network formation is stabilized by the positive 
tensions. 

When a positive tension 7 is imposed for the coex¬ 
isting states of network and tubules as shown in the 
second snapshot of Fig. □Id), the tubules continue to 
grow at irfodAsT = 40. However, the tubules shrink 
at 7 r 2 od /fcBT = 70. Thus, at the critical tension 
7 r rod/^BF — 50, the tubule elongation force / t b by the 
rod assembly is balanced with the expansion of the pro¬ 
jected membrane area by the surface tension as ftb = 
27ri?tb7- This tension is higher than that required to 
suppress the tubule protrusion from the flat membranes 
(iWod/^B ^ 1 — 30), because a nucleation barrier exists for 
the protrusion. 




3 



FIG. 1: Tubulation dynamics from flat membrane for high rod density, </> ro d = 0.4 (A^od = 1,024). (a) Protein rod with 
spontaneous rod and side curvatures, C ro d and C s id e , respectively. The protein rod is displayed as a chain of spheres, the halves 
of which are coloured red and yellow. The orientation vector lies along the line of the yellow to red hemispheres. The light 
blue spheres represent membrane particles, (b) Dynamic phase diagram of tubulation from a tensionless flat membrane. The 
red circles represent percolated network formation before tubulation. The green diamonds indicate that the tubules are formed 
from isolated clusters. The blue squares represent partial network formation, and the dashed lines are guides for the eye. (c,d) 
Sequential snapshots of tubulation from tensionless flat membrane for Csidefrod = ( c ) 1 and (d) —1 at Crod^rod = 4. (e-g) Time 
evolution of mean cluster size (N c i) and mean cluster height ((z^)) 1 / 2 for Gside^rod = (e) 1, (f) 0, and (g) —1 at Crodn-od = 4. 
Error bars calculated from eight independent runs are displayed at several data points. 
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FIG. 2: Snapshot of tubulation from tensionless flat mem¬ 
brane at C ro dr ro d = 2.5, Cside^rod = -1, and (f ) IO d = 0.4. 


Geometrical Analysis 

To clarify the effects of C s id e on the network forma¬ 
tion, the difference between the free energy of a hexago¬ 
nal array of the rod assembly and that of a striped array 
is estimated using a simple geometric model (see Fig. [4]). 
Here, a percolated network is modelled as a hexagonal ar¬ 
ray with side length [Fig. SJa)], while an unbranched 
rod assembly is modelled as a striped array [Fig. 01(c)] . 
The rod assemblies have rectangular shapes with widths 
equal to the rod length r ro( j. Our analysis shows that the 
hexagonal network can have lower energy for C s ;d e < 0. 
This explains why the membranes are trapped in the 
branched network as a local free-energy minimum. 

In the striped array, the rod assemblies are aligned in 
parallel with intervals of (1 — (/w^rod/^rod- The rod 
assemblies are curved upwards and the other regions are 
curved downwards, as shown in Fig. 01(c). To maintain 
the continuity of the normal vector of the membrane, the 
curvatures have the relation C 2 = </> r odCi/(l —0rod)- The 
curvature energy F at of the striped array per the area r 2 od 
is given by 

-f L = ^(Cr-CV odf^rod ( 1 ) 

r rod 2 

+ ~C s 2 ide<^rod + - C\ (1 — (j) I0 d), 

where K rl and k i2 are the bending rigidities of the rod 
assembly parallel and perpendicular to the rod axis, re¬ 
spectively. In our simulation, k = 15kBT, K r i = AOk^T, 
and k T 2 = k [42j]. 

In the hexagonal array case, the membrane surface 
is divided into the following three regions: Region I: 
A rectangular rod assembly with length L h — v^ryod/S 
and width r ro d! Region II: A triangular membrane with 
side length r ro d at the vertex; and Region III: A hexag¬ 
onal membrane with side length Lh — %/3r ro d/3. On 
a flat membrane, the areas of these regions are given 
by A\ = (L h - v / 3r ro d/3)r ro d, An = V3r* od /A, and 
Am = 3-\/3(Th — v / 37 'rod/ 3 ) 2 / 2 , respectively. To sim- 



FIG. 3: Membrane tubulation from flat membrane for low 
<Aod = 0.1 (iV ro d = 256), at C ro d»~rod = 4. (a) Sequential 
snapshots at Cside = 0 and 7 = 0 . (b) Time evolution of 
(Ad) and ((zj?i))h 2 . The solid lines represent the data for 
Cside^rod = —1, 0, and 1 at 7 = 0. The gray and black 
dashed lines represent the data for -yr^ od /kBT = 25 and 37.5 
at Cside = 0, respectively. Error bars calculated from eight 
independent runs are displayed at several data points. 


plify the calculation, it is assumed that each region 
has constant curvatures and its area is independent of 
these curvatures. Region I has curvatures Ci along the 
rod axis and —C 3 perpendicular to the rod axis [along 
the dashed line in Fig. 01(a)], such that it has a saddle 
shape. Regions II and III are triangular and hexag¬ 
onal spherical caps with radii l/Ci and I/C 4 , respec¬ 
tively. To maintain the continuity of the normal vector 
of the membrane, C 3 = (tt — 20 2 )/(L\ v — V3r IO< y/3) and 
C 4 = Ci7\ od /(V3L h - r ro d)- The angle 6 2 is given by 
cos 2 (0 2 ) = (1/ cos 2 ($i/2) — l)/3, where 9y = Cir ro d, 
since the vertices of three rod assemblies make contact 
with each other on the spherical cap and maintain three¬ 
fold rotational symmetry. At 0 2 = 0, Ci has a maximum 
value of 27r/3r ro d. As the area fraction of region I is 
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FIG. 4: Energy analysis of network structure using simple 
geometric model, (a-c) Schematic representation of geometric 
model, (a) Top view of hexagonal array at Ci = 0. (b) Side 
view of hexagonal array along the dashed line in (a), (c) 

Bird’s-eye and front views of striped array, (d) Free-energy 
difference A F between hexagonal and striped arrays of rod 
assembly. The solid lines represent the data for Csiderrod = 
— 1, 0, and 1 at (f> IO A = 0.4. The dashed line represents the 
data for Cside^rod = — 1 and cf > IO d = 0.5. 


cj >rod = 2v / 3v4i/3L^, Lh is obtained as 


Lh 


2r rod 

V3(l - v / T^2Q ' 


( 2 ) 


The curvature energy Lh ex of the hexagonal array per 

r rod is S iven by 

%= ^-(Ci-C r od )V rod (3) 

Cod 2 

+ + Csi d e) 2 </>ro d 

+ kCi(1 — </> ro d ~ y/l — 2 (f> I0 d) 

+ kC|( 1 — 0 rod + \] 1 — 2 (f> ro d)- 



FIG. 5: Snapshots of vesicles at (p IOd = 0.3 and N io a = 
288. (a) Spherical shape in thermal equilibrium at C ro d = 
0 and Gside = 0. (b) Elongated discoidal shape in thermal 
equilibrium at C ro dr rod = 4 and C s ider ro d ■ — 1- (c) Three 
metastable shapes at C ro dr rod = —4 and C s id e r r od = 1- The 
membrane particles are displayed as small transparent spheres 
for clarity. 


In our simulation, networks with a wider range of C s i de 
are formed for smaller values of C' rod [see Fig.QJb)]. This 
dependence can be explained by an effective increase in 
the area fraction </> rod . The rod assembly region contains 
more membrane particles with decreasing C'rod (see light 
blue particles in rod network and tubules in Fig. EH). 
Thus, the area of region I increases. At </> rod = 0.5, the 
minimum of A F is twice that at (j) ro d = 0.4, although the 
area fraction is only 25% larger [see the dashed line in 
Fig. 01(d)] . Thus, the effective increase in the rod region 
enhances the network formation. 

As <^ rod decreases, Lh of the hexagonal array increases 
[see equation ©]• However, it is difficult for the hexago¬ 
nal network with long Lh to form spontaneously, since the 
tubule formation begins before the rod assembly reaches 
Lh- The simulation results for </> rod =0.1 indicate tubu- 
lation from clusters with N c \ > 40 at C ro dC-od = 4. 
Thus, the formation threshold of the percolated network 
is Lh — 5r rod) where ^ rod — 0.2. In the simulation at 
4>Tod = 0.2, the clusters are typically percolated only in 
one direction only at C rod r r od = 4 and C s ide^rod = — 1, 
which supports this estimation of the critical density. 


Shape Transformation of Small Vesicles 


The energy difference A F = Lh ex — F st is shown in 
Fig.[4](d). The first terms in equations m and m cancel 
since both arrays have the same curvature C\ along the 
rods. Hence, A F is independent of C'rod and k t \ . In the 
hexagonal array, the curvature energy of the rod assem¬ 
bly is reduced for C s i de < 0, and Lhex is smaller than F st 
for small C\. Thus, the branched network can be stabi¬ 
lized by a negative Cside during the formation of the rod 
assembly. 


Next, we investigate the tubulation from a vesicle 
of radius R ves = 3.07r rod at (j) TO d = 0.3 (see Fig. [5]) 
and demonstrate that the original membrane curvature 
C ves = 1 / Rves = 0.33/r ro d changes the tubulation dy¬ 
namics. For the positive curvature C ro d?~rod = 4, no tubu¬ 
lation is obtained for C s ide^rod = — 1, 0, and 1. Instead, 
the vesicle deforms to an elliptic disk and the rods sur¬ 
round the disk rim [see Fig.[5](b)]. A discoidal bud is often 
transiently formed, but the rearrangement of the rod as- 












6 


semblies results in disk formation even for Cside^rod = — 1 
(see Supplemental Movie 5). Thus, outward tubulation 
and network formation are suppressed in small vesicles. 

In contrast, for the negative curvature C ro d^rod = — 4, 
tubulation into the inside of the vesicle is obtained (see 
Fig. Etc) and Supplemental Movie 6 ). A percolated net¬ 
work with CsideFrod < 0 has a significantly longer lifetime 
than in the flat membrane case [see the middle and right 
snapshots in Fig. Etc)] . The coexistence of tubules and a 
ring is also obtained [see the left snapshot in Fig. Etc)]. 
The ring stabilizes an outward bud. For C ro d < 0, the 
rods bend the membrane towards the interior (opposite 
to the original membrane curve), such that the rods lo¬ 
cally form a saddle shape in which the two principal cur¬ 
vatures have opposite signs. The network and ring struc¬ 
tures are stabilized by the positive (opposite) C s id e > 0 , 
but not by C s ide < 0 . 


Discussion 

We have revealed that, in addition to the spontaneous 
curvature along the protein rods C ro d, the perpendicu¬ 
lar spontaneous curvature C s ide significantly influences 
the protrusion of membrane tubules. The percolated- 
network structure of the rod assembly has a long lifetime 
for C s ide < 0 , because the saddle membrane shape at 
branches of the rod network is stabilized by the opposite 
curvature of C s jd e with respect to C ro d- Thus, the net¬ 
work formation decelerates the tubulation significantly, 
despite having a minor effect on the equilibrium property. 
Both positive surface tensions and membranes originally 
bending in the same direction as C s jd e < 0 can stabilize 
the network structure. Our findings provide new insights 
into the regulation of biomembrane shapes by curvature- 
inducing proteins. 

Here, we employ Langevin dynamics, in which hydro- 
dynamic interactions are neglected. Since the static sta¬ 
bility of the network branch is the key factor, we do not 
expect the obtained C s ide dependence to be qualitatively 
changed by the hydrodynamic interactions. However, the 
network formation condition may be modified. The diffu¬ 
sion coefficient of the proteins on the membrane depends 
on the protein size, and fast protein diffusion compared 
to the membrane deformation speed likely enhances the 
network formation. 

The F-BAR domain of Pascin is considered to have a 
nonzero side curvature, since it has an S-shape on the 
membrane surface in the addition to the curvature per¬ 
pendicular to the membrane Jlol |. Pascin induces mem¬ 
brane tubes in a wide range of diameter In the 

present simulations, we did not obtain such a behaviour. 
It may require a larger side curvature or attractive inter¬ 
actions between the rods. The rod assembly with large 
side curvatures is an interesting problem for further stud¬ 
ies. 


An assembly of F-BAR proteins, Cdcl5, has been ob¬ 
served along the contractile ring of cell division 44]. 
Their adsorption to the inner leaflet of the plasma mem¬ 
brane is considered to yield a ring structure similar to 
that shown in Fig.[5](c). Our study suggests that the side 
curvature may play an important role in the formation 
of neck-like structures during cell division and membrane 
budding in endo/exocytosis. 


Methods 


We employ one of the meshless membrane models 
f4l| . in which a fluid membrane is represented by a 
self-assembled one-layer sheet of membrane particles. A 
membrane particle has an excluded volume with diameter 
a and an orientational degree of freedom. The solvent is 
implicitly accounted for by an effective potential between 
the membrane particles. The mechanical properties of 
the fluid membrane can be varied over a wide range. The 
details of the meshless membrane model and protein rods 
are described in Ref. 41 and Ref. HH respectively. In this 
study, we employ the parameter set used in Ref. HI for a 
membrane with isotropic spontaneous curvature Co = 0. 
The membrane has mechanical properties typical for lipid 
membranes: Bending rigidity n/k^T = 15 ± 1, tension¬ 
less membrane area per particle ao/a 2 = 1.2778±0.0002, 
area compression modulus K^a 2 /ksT = 83.1 ± 0.4, and 
edge line tension Ta/k^T = 5.73 ± 0.04. 

A BAR protein is modelled as a curved rod consisting 
of a chain of N sg membrane particles with r ro( j = 10er 
and N sg = 10. The rod has anisotropic spontaneous 
curvature C ro d along its length and spontaneous curva¬ 
ture C s ide perpendicular to its length. When two protein 
rods come into contact, C s id e is applied between the rods. 
When the rod is surrounded by membrane particles, the 
spontaneous curvature C s id e /2 is applied between the rod 
and neighbouring membrane particles. A molecular dy¬ 
namics with a Langevin thermostat is employed 41], |45|. 
The simulation results are displayed with a time unit of 
r = r 2 od /D , where D is the diffusion coefficient of the 
membrane particles in the tensionless membranes. We 
use total particle numbers N = 25, 600 and 9, 600 for flat 
membranes and vesicles, respectively. The rod density is 
defined as </> ro d = N Iod N sg /N. 

The mean square cluster height (z 2 ^ is calculated 
as follows. A rod is considered to belong to a clus¬ 
ter when the distance between the centres of mass of 
the rod and one of the rods in the cluster is less than 
r ro d/2. The height variance of each cluster is calculated 
as z i, c \ = J2k( z k ~ z itCm ) 2 /Ni^iN sg , where N itc i is the 
number of rods belonging to the z-th cluster, 2 j jCm is the 
z component of the centre of mass of the cluster, and 
the summation is taken over all rod segments in the clus¬ 
ter. Finally, z 2 1 is calculated as the average of zf cl for all 
clusters. 
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